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ABSTRACT 

Circuit  theory  is  applied  in  the  determination  of  the  time-varying  axial  magnetic  field  produced 
by  the  interaction  of  a  vibrating  aluminum  enclosure  and  the  earth's  magnetic  field.  The  axial 
magnetic  field  that  exists  within  the  enclosure  is  received  by  a  loop  antenna.  A  steady-state  time- 
harmonic  analysis  is  applied.  Theoretical  results  are  compared  with  measured  data  at  100  Hz. 
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1.  INTRODUCTION 

Consider  an  aluminum  enclosure  of  length  L  =  76  in.  with  a  triangular  cross  section  of  width 
W  =  33  in.  and  height  H  =  25  in.  as  shown  in  Fig.  1.  The  walls  of  the  enclosure  are  4  in.  thick 
with  the  exception  of  the  base  which  is  5  in.  thick.  The  enclosure  is  sealed  at  both  ends  and 
provides  sufficient  shielding  of  atmospheric  noise  for  frequencies  of  50  -  100  Hz.  Within  the 
enclosure  lies  a  loop  antenna  with  a  magnetic  core.  The  loop  antenna  is  in  the  form  of  a  solenoid 
consisting  of  30,800  turns  of  26  gauge  copper  wire  that  are  wound  in  6  separate  sections  with  22 
layers  per  section.  The  magnetic  core  is  a  circular  cylinder  of  4  ft  length  and  2  in.  diameter. 

The  enclosure  is  made  to  vibrate  while  the  loop  antenna  is  held  stationary  within  a  vibration 
isolation  system.  The  interaction  of  the  moving  enclosure  with  the  earth's  dc  magnetic  field  is 
equivalent  to  a  time-varying  magnetic  flux  within  the  enclosure,  inducing  a  current  in  the 
aluminum.  The  component  of  current  flowing  in  the  cross-sectional  plane  of  the  enclosure 
produces  a  time-varying  axial  magnetic  field  within  the  enclosure  that  is  received  by  the  loop 
antenna.  The  frequency  of  the  measured  magnetic  field  is  correlated  with  the  vibration  frequency 
of  the  enclosure. 

The  purpose  of  this  study  is  to  determine  the  time-varying  axial  magnetic  field  received  by  the 
loop  antenna  due  to  the  vibratory  motion  of  the  aluminum  enclosure.  A  steady-state  time-harmonic 
analysis  will  be  applied. 


2.  ENCLOSURE  MODEL 

Fig.  2  illustrates  three  types  of  vibratory  motion  of  the  enclosure:  (a)  translational  motion  in  the 
radial  direction,  (b)  translational  motion  in  the  axial  direction,  and  (c)  rotational  motion.  The  actual 
vibratory  motion  consists  of  a  combination  of  these  types  of  motion.  The  enclosure  is  made  of 
aluminum,  a  good  conductor  (a  =  3.77  x  107  S/m),  and  can  be  thought  of  as  a  fluid  of  free 
electrons.  As  the  enclosure  vibrates,  the  free  electrons  in  the  metal  interact  with  the  earth's  dc 
magnetic  field  resulting  in  a  Lorentz  force  experienced  by  each  electron.  The  net  effect  of  these 
forces  is  to  produce  an  electric  current  in  the  aluminum.  Because  of  the  vibrating  motion  of  the 
enclosure,  this  current  will  be  time  varying. 

A  current  circulating  in  the  cross-sectional  plane  of  the  enclosure  is  necessary  to  produce  an 
axially  directed  magnetic  field  that  can  be  received  by  the  loop  antenna.  Therefore,  the  only  motion 
of  interest  in  this  problem  is  one  which  results  in  a  circumferentially  directed  current.  To  better 
understand  which  types  of  vibratory  motion  will  produce  a  circumferential  current  in  the 
aluminum,  let  us  approximate  the  geometry  of  the  enclosure  by  a  tubular  cylinder  as  in  Fig.  2.  For 
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the  translational  motions  in  Fig.  2  (a)-(b),  diametrically  opposite  electrons  in  the  cylinder  will 
experience  the  same  Lorentz  force,  resulting  in  a  current  distribution  with  no  net  circumferential 
flow  regardless  of  the  orientation  of  the  earth's  magnetic  field.  Consequently,  translational  motion 
is  of  no  interest  in  this  problem.  In  contrast,  rotational  motion  of  the  enclosure  results  in  a  time- 
varying  axial  magnetic  flux  through  the  shield,  producing  an  induced  electromotive  force  (emf)  in 
the  enclosure.  This  induced  emf  produces  a  time-varying,  circumferentially  directed  current 
distribution  in  the  aluminum. 

To  determine  the  induced  emf  in  the  enclosure,  consider  the  rotational  vibration  of  the 
enclosure  about  its  center  in  the  x-z  plane  with  coordinate  system  shown  in  Fig.  3,  where  the  z  axis 
is  coincident  with  the  axis  of  the  enclosure.  The  primed  coordinates  refer  to  the  rest  frame  of  the 
enclosure  while  the  unprimed  coordinates  refer  to  the  earth's  reference  frame.  The  angle  0  that  the 
primed  coordinate  axes  make  with  respect  to  the  unprimed  axes  is  given  by 

0(t)  =  0O  sin  cot  (1) 

where  co  =  27cf  is  the  angular  frequency  of  vibration  and  0O  is  the  amplitude  of  vibration.  Note  that 
for  vibrations  of  interest,  I0ol  «  1 . 

The  earth's  dc  magnetic  field  (flux  density)  is  given  in  each  reference  frame  as 

B  =  i  Bxo  +  j  By0  +  k  Bzo  (2a) 

B'  =  i'  Bx  +  j'  By  +  k'  Bz  (2b) 


where  i,  j,  and  k  are  the  unit  vectors  along  the  x,  y,  and  z  axes,  respectively,  while  i',  j',  and  k' 
are  the  unit  vectors  along  the  x',  y',  and  z'  axes,  respectively.  Since  the  vibration  is  confined  to  the 
x-z  plane,  j  =  j'.  The  components  of  the  earth's  field  in  the  enclosure  frame  are  given  in  terms  of 
the  unprimed  components  as  follows: 


$ 

Bx  =  Bxo  cos  0  +  Bzo  sin  0 

(3a) 

By  =  Byo 

(3b) 

Bz  =  -Bxo  sin  0  +  Bzo  cos  0 

(3c) 

Note  that  the  field  components  in  the  primed  frame  in  eqs.  (3)  are  time  varying  since  0  is  a  function 
of  time.  The  axial  magnetic  flux  passing  through  the  enclosure  is 


0(t)  = 


A  sin  0  +  Bzo  A  cos  0 


(4) 
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where  A  is  the  cross-sectional  area  within  the  enclosure.  It  should  be  observed  that  the  y 
component  of  the  earth's  magnetic  field  produces  no  contribution  to  the  axial  magnetic  flux.  The 
rate  of  change  of  flux  through  the  enclosure  is 

O'(t)  =  -BxoA(o0o  coscot  cos0(t)  -  BzoAcoOo  coscot  sin0(t) 

2 

=  -BxoAco0o  coscot  -  BzoA^0o  sin2cot 


=  -  Re  (Viei®1  +  V2ei2a)t) 

(5) 

where 

V  i  =  BxoAa>0o 

(6a) 

V2  =  -jBzoA^0o 

(6b) 

In  (5)  the  approximation  is  valid  provided  I0ol  «  1.  Vj  and  V2  are  induced  emf  phasors 
corresponding  to  frequencies  co  and  2co,  respectively.  Note  that  V  j  is  generally  dominant  since 
I0ol  «  1.  Consequently,  the  x  component  of  the  earth's  magnetic  field  produces  the  most 
significant  ac  current  in  the  enclosure  at  the  same  frequency  as  the  vibration.  The  above  expression 
for  Vj  will  be  used  in  the  next  section  to  determine  the  current  distribution  in  the  enclosure. 

3.  CIRCUIT  ANALYSIS 

The  vibrating  enclosure  can  be  viewed  as  a  single-turn  loop  with  a  circuit  representation  as  a 
lossy  inductor  in  the  frequency  range  of  interest  (50  -  100  Hz).  Fig.  4  illustrates  a  simple  circuit 
model  showing  the  interaction  of  the  vibrating  enclosure  and  the  enclosed  antenna.  Note  that  the 
enclosure  and  the  loop  antenna  are  inductively  coupled  through  the  mutual  inductance  M.  The 
circuit  equations  for  Fig.  4  are  as  follows: 

v(t)  =  Rii1+(Li  +  Le)^-  +  M^  (7a) 

0  =  M^  +  (Ra  +  RL)i2  +  La^  (7b) 

where  R;  and  L;  are  the  internal  resistance  and  inductance  of  the  enclosure,  respectively,  while  Le 
is  the  external  inductance  of  the  enclosure.  Ra  and  La  represent  the  resistance  and  inductance, 
respectively,  of  the  loop  antenna  while  Rl  is  the  load  resistance.  The  emf  voltage  induced  in  the 
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enclosure  is  given  by  v(t)  while  ij(t)  is  the  total  current  flowing  circumferentially  in  the  enclosure. 
The  current  i2(t)  flows  through  the  loop  antenna.  The  dot  convention  in  Fig.  4  refers  to  the 
calculation  of  the  mutual  inductance  between  the  enclosure  and  loop  antenna  [1]. 

Next  we  will  consider  the  steady-state  frequency  response  of  the  circuit  in  Fig.  4  given  an 
excitation  frequency  of  co  rad/s.  If  the  currents  and  voltage  in  Fig.  4  are  replaced  by  phasor 
quantities,  the  circuit  equations  (7)  become 

V  =  [Rj  +  jco  (Li  +  Le)]  Ii  +  jcoMIz  (8a) 


0  —  jcoMIi  +  (Ra  +  Rl  +  jcoLa)  L 

The  solution  of  (8)  for  the  phasor  current  It  is 

j  _  Y _ Ra  +  Rl  +  jcoLa _ 

co2M2  +  [  R}  +  jco  (Li  +  Le)]  [  Ra  +  RL  +  jcoLa] 


(8b) 


(9) 


Given  a  steady-state  sinusoidal  emf  voltage  V,  the  total  phasor  current  It  flowing  circumferentially 
in  the  enclosure  is  given  in  (9).  In  this  analysis  the  current  induced  in  the  enclosure  is 
approximated  by  a  uniform  solenoidal  current  sheet  with  any  end  effects  ignored.  The  end  effects 
can  be  safely  ignored  provided  the  length  of  the  enclosure  is  much  greater  than  the  maximum 
cross-sectional  dimension. 

To  simplify  the  determination  of  the  circuit  parameters  of  the  enclosure,  consider  the 
geometrical  approximation  of  the  actual  enclosure  by  a  cylindrical  tube  of  length  L,  inner  radius  b, 
and  thickness  t  as  shown  in  Fig.  5.  The  radius  b  is  determined  by  equating  the  inside  cross- 
sectional  areas  of  the  enclosure  and  the  tubular  model,  yielding 

-i IW 

where  H'  and  W'  are  the  inside  height  and  width  of  the  enclosure,  respectively.  From  the 
enclosure  dimensions  in  Fig.  1,  we  find  that  H'  =  13.7  in.  and  W'  =  22.6  in.  Therefore,  (10) 
yields  b  =  7.03  in.  The  axis  of  the  loop  antenna  is  assumed  to  be  coincident  with  the  axis  of  the 
tubular  enclosure  in  Fig.  5. 


3.1  CURRENT  DISTRIBUTION 

The  current  density  in  the  tube  is  given  approximately  in  cylindrical  coordinates  by  [2] 
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J  =  a<p  J0  e*(  1+j)(p-b)/5  ,  b  <  p  ^  b  + 1 


(11) 


where  8,  the  skin  depth  in  aluminum,  is 

8  =  -=J= 

Yrcf |i0c 


(12) 


In  (1 1),  by  is  the  unit  vector  in  the  azimuth  direction,  p  is  the  radial  coordinate,  and  JQ  is  the 
current  density  along  the  inner  surface  of  the  tube.  In  (12),  f  is  the  vibration  frequency  of  the 
enclosure  and  (Iq  is  the  permeability  of  free  space.  Note  that  the  following  approximations  have 
been  made  in  (11): 

(a)  The  current  density  is  the  same  as  that  induced  in  a  plane  conductor.  This  is  valid  for  the 
tubular  enclosure  provided  b/8  >7  [2]. 

(b)  The  current  distribution  is  independent  of  <p.  This  is  valid  since  2nb  «  X,0,  where  XQ  is  the 
wavelength  in  free  space. 

(c)  J  is  independent  of  z  as  end  effects  have  been  ignored.  This  is  valid  since  L/b  »  1. 

The  total  current  flowing  in  the  enclosure  is  given  by 


Ii 


Jo 


e-(i+j)(p-b)/S  dp  dz 


=  J0L 


^L_[l  .  e-(l+j)t/5] 


(13) 


If  t  >  48,  (13)  reduces  to 


1 1  =  J0L 


_ 8 _ 

(1  +  j) 


(14) 


From  (11)  and  (14),  the  current  density  J  can  be  written  in  terms  of  the  total  current  Ij  as 

j  =  a  klLtfi  e-(H)(P-b)/S  s  b  <  p  <  b  +  t 

L  8 


(15) 


To  simplify  the  analysis,  the  current  distribution  may  be  approximated  as  a  solenoidal  current  sheet 
and  given  by 

J  =  a<pi-8(p-beff)  (16) 

where  8  is  the  Dirac  delta  function  and  beff  is  the  effective  radius  defined  by 
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beff  = 


J("b+t 


p  e-(P-b)^>  dp 

b  =  -t5  e~tj/s  +  (5+b)5(l-e~tj/s) 

5(l-e-^) 


Jpb+t 

„ 


e-(p-b)/5  dp 


If  t  >  45,  (17)  reduces  to 


(17) 


beff  =  b  +  5  (18) 

Thus  the  induced  current  in  the  tubular  enclosure  is  approximated  by  a  circumferentially  directed 
surface  current  Ij/L  that  is  uniformly  distributed  over  its  length  with  an  effective  radius  of  b  +  5. 
This  approximation  will  simplify  the  determination  of  the  external  inductance  Le  of  the  enclosure 
and  the  mutual  inductance  M  between  the  loop  antenna  and  the  enclosure. 


3.2  DETERMINATION  OF  CIRCUIT  PARAMETERS 

The  internal  resistance  R;  and  inductance  L;  of  the  enclosure  can  be  determined  from  the 
expression  for  the  internal  impedance  of  a  flat  plane  conductor  of  length  2jtb  and  are  given  by  [2] 

Ri  =  coLi  =  2£lL  (19) 

o5L 


This  formula  is  valid  provided  L  »  b  and  t  >  35. 

The  external  inductance  Le  of  the  tube  is  determined  by  approximating  it  as  a  solenoidal  current 
sheet  of  radius  beff  as  expressed  by  (16).  In  this  case  the  inductance  is  given  by  [3] 

<20) 

where  K  is  Nagaoka's  constant,  a  function  of  the  radius  and  length  of  the  sheet.  K  is  given  by  a 
complicated  relation  involving  elliptic  integrals  and  is  plotted  in  [4],  The  expression  (20)  assumes 
that  the  current  sheet  encloses  a  region  that  is  filled  with  air  (permeability  =  (ij.  However,  since  a 
loop  antenna  with  a  magnetic  core  of  permeability  (i  lies  within  the  enclosure,  the  effect  of  this 
medium  must  be  accounted  for  in  the  determination  of  Le.  The  magnetic  core  may  be  accounted  for 
by  replacing  (i0  in  (20)  by  an  effective  permeability  (ieff.  Thus  (20)  is  replaced  by 
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T  _|ieff7Cb^ffK 
_  " 


(21) 


where 


(22) 


In  (22),  a  is  the  radius  and  (ij  is  the  relative  permeability  of  the  core. 

The  mutual  inductance  M  between  the  loop  antenna  and  the  enclosure  can  be  obtained  by 
considering  the  coupling  between  two  concentric,  coaxial,  solenoidal  current  sheets  as  described  in 
Fig.  6.  If  the  loop  antenna  is  approximated  as  a  single  layer  solenoid,  the  mutual  inductance  is 
approximately  given  by  [5] 


M=  1.57  ^eff^  — 1 


1  +&4(3  _4)1 


8X 


4  V 


(23) 


where 

x=Va'+^  (24) 

In  (23),  Lj  is  the  axial  length,  aj  is  the  radius,  and  Nj  is  the  total  number  of  turns  on  the  outer 
solenoid  (enclosure)  while  L2,  a2,  and  N2  are  the  corresponding  quantities  for  the  inner  solenoid 
(loop  antenna).  The  effect  of  the  magnetic  core  of  the  loop  antenna  is  accounted  for  through  the 
effective  permeability  as  defined  in  (22).  Expression  (23)  is  an  extension  of  a  formula  derived 
for  the  case  of  an  air  core.  Additional  formulas  similar  to  (23)  also  appear  in  the  literature  [6  -  7] 
for  the  case  of  an  air  core.  In  the  application  of  (23)  to  this  problem,  aj  is  the  effective  radius  beff 
of  the  enclosure  as  defined  in  (18)  while  a2  is  approximated  as  the  average  of  the  inner  and  outer 
radii  of  the  loop  antenna. 


4.  DETERMINATION  OF  THE  AXIAL  MAGNETIC  FIELD 

Consider  a  current  distribution  J  in  an  infinite  homogeneous  medium  as  described  in  Fig.  7. 
From  the  Biot-Savart  law  [8],  the  total  magnetic  field  at  a  point  r  is  given  by 
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H(r)=-L 

4n 


J(r')  x  (r  -  r') 


r  -  r 


/1 3 


dv' 


(25) 


where  the  primes  denote  integration  variables  while  r  and  r'  denote  position  vectors  corresponding 
to  the  observation  and  source  points,  respectively.  The  geometry  for  application  of  the  Biot-Savart 
law  is  given  in  Fig.  8.  The  application  of  the  Biot-Savart  law  (25)  for  the  uniform  solenoidal 
current  sheet  given  in  (16)  yields  the  following  surface  integral: 


H(r)=4^L 


f  , 

a,x(r- r') 


r  -  r 


A  3 


■dA' 


(26) 


where  the  integration  is  over  the  current  sheet  of  length  L  and  radius  beff.  The  axial  component  of 
magnetic  field  along  the  axis  of  the  tube  at  (0,0,z)  is  determined  through  evaluation  of  (26)  in 
cylindrical  coordinates,  yielding 


r  l 


Hz(0,0,z) 


h_ 

2L 


(27) 


From  [9],  No.  165, 


_dx_ 


(x2  +  a2)312  a2  Vx2  +  a2 


(28) 


where  a  is  a  constant.  The  application  of  (28)  to  (27)  yields 


Hz(0,0,z)  = 


L-  z 


V  Z2  +  bgff  V(L  -  zf  +  bgff 


(29) 


The  magnetic  field  at  the  center  of  the  tube  (0,0,172)  is 


Hz(0,0,L/2)  = 


Ii 


V  L2  +  4bl 


eff 


(30) 


Note  that  the  effect  of  the  magnetic  core  has  been  ignored  in  (25)  -  (27).  As  a  first  approximation, 
the  magnetic  field  at  the  center  of  the  tube  will  not  change  appreciably  from  (30)  when  the  magnetic 


8 


TM  931044 


core  is  present  since  Hz  is  continuous  across  the  longitudinal  surface  of  the  core. 


5.  NUMERICAL  RESULTS 

In  this  section,  we  will  determine  the  axial  magnetic  field  produced  by  the  enclosure  at  a 
vibration  frequency  of  100  Hz.  In  the  vicinity  of  the  measurement  site  at  Fishers  Island,  NY,  the 
earth's  magnetic  flux  density  is  B  =  0.5  Gauss  =  5  x  10'5  Wb/m2  at  an  angle  of  approximately  70° 
with  respect  to  the  axis  of  the  enclosure  as  shown  in  Fig.  9.  It  is  assumed  that  the  earth's  magnetic 
induction  vector  B  lies  entirely  in  the  plane  of  vibration  of  the  enclosure.  Consequendy, 

Bxo  =  B0  cos  70°  =  1.71  x  lO  5  Wb/m2. 

This  assumption  will  result  in  a  maximum  ac  magnetic  field  produced  by  the  enclosure.  From 
(6a),  the  induced  emf  voltage  (per  radian  of  peak  vibration  0O)  in  the  enclosure  is 

V  =  Vj  =  Bxojtb2co  =  1.08  x  10-3  Volts/rad. 

Note  that  this  voltage  is  at  100  Hz. 

At  100  Hz,  the  skin  depth  in  aluminum  is  (from  (12))  8  =  0.323  in.  Since  the  thickness  (t  =  4 
in.)  of  the  aluminum  is  greater  than  4  skin  depths  (t/8  =  12.4),  the  approximation  (18)  is  valid, 
yielding  an  effective  enclosure  radius  beff  =  b  +  8  =  7.35  in.  The  conditions  assumed  in  section 
3.1  are  satisfied  (b/S  =  21.8  >  7,  27tbA0  =  3.74  x  10'7  «  1,  L/b  =  10.8  »  1,  t/8  =  12.4  >  4), 
thereby  validating  the  approximation  of  the  current  density  given  by  (15). 

The  loop  antenna  has  a  magnetic  core  with  effective  permeability  m  =  230  and  radius  a  =  1  in. 
Therefore,  the  effective  permeability  of  the  medium  interior  to  the  enclosure  is  peff  =  5.24  p0  as 
obtained  from  (22).  With  L  =  76  in.  and  beff  =  7.35  in.,  the  length  to  effective  diameter  ratio  of  the 
enclosure  is  L/2beff  =  5.17.  Consequently,  from  [4],  Nagaoka's  constant  for  the  enclosure  is  K  = 
0.92,  resulting  in  an  external  inductance  Le  =  0.344  pH  as  obtained  from  (21).  From  (19),  the 
internal  resistance  and  reactance  of  the  enclosure  are  Rj  =  coL;  =  1.88  pf2.  Note  that  the  external 
reactance  coLe  (=216  pf2)  is  much  greater  than  the  internal  resistance  and  reactance. 

In  the  determination  of  the  mutual  inductance  M  between  the  enclosure  and  the  loop  antenna, 
the  parameters  of  the  enclosure  are  Nj  =  1  turn,  aj  =  beff  =  7.35  in.,  and  Lj  =  L  =  76  in.  The 
parameters  for  the  loop  antenna  are  N2  =  30,800  turns,  a2  =  1.30  in.,  and  L2  =  24.25  in.  Note 
that  a2  is  the  average  of  the  inner  and  outer  radii,  1.125  in.  and  1.475  in.,  respectively,  of  the  loop 
antenna.  The  substitution  of  these  parameters  into  (23)  yields  M  =  353  pH. 

The  circuit  parameters  of  the  loop  antenna  as  shown  in  Fig.  4  are  Ra  =  845  £2,  La  =  592  H,  and 
Rl  =  15  MQ.  These  values  were  determined  experimentally.  With  the  above  circuit  parameters,  we 
obtain  the  following  impedances  at  100  Hz:  coM  =  2.22  m£2,  IRj  +  jco(L;  +  Le)l  =  coLe  =  0.216  m£2, 
IRa  +  Rl  +  jcoLJ  =  Rl  =  15  M£2.  Therefore, 
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(coM)2  «  |R,  +  jco(Li  +  Le)|  |Ra  +  Rl  +  jcoLa| 

Thus (9) reduces  to 

I  ~  V _ =  _Y_  /3n 

1  R,  +  jco(L,  +  Le)  jcoLe  V 

The  above  expression  indicates  that  the  enclosure  is  effectively  decoupled  from  the  loop  antenna. 
From  (31),  the  magnitude  of  the  total  current  per  radian  of  peak  vibration  0O  in  the  enclosure  is 
II] I  =  4.95  A/rad.  Therefore,  from  (30),  the  magnetic  field  strength  per  radian  of  peak  vibration  0O 
at  the  center  of  the  loop  antenna  is  IHZ(0,0,L/2)I  =  2.52  A/m-rad.  From  experiment,  the  peak 
vibration  angle  is  0O  =  1.40  x  10'8  rad.  Thus  the  magnetic  field  strength  is 

IHz(0, 0,L/2)I  =  3.53  x  108  A/m  =  -149  dBA/m. 

In  comparison,  the  field  strength  measured  by  the  loop  antenna  was  -140  dBA/m,  9  dBA/m  higher 
than  the  computed  result.  The  difference  between  the  computed  and  measured  results  may  be 
attributed  to  the  various  approximations  made  in  the  model.  In  spite  of  this  difference,  the  analysis 
shows  that  the  vibrating  enclosure  produces  a  significant  ac  magnetic  field  at  100  Hz. 


6.  SUMMARY  AND  CONCLUSIONS 

A  steady-state  time-harmonic  analysis  has  been  applied  in  the  determination  of  the  time-varying 
axial  magnetic  field  produced  by  a  vibrating  enclosure  in  the  presence  of  the  earth's  magnetic  field. 
The  analysis  was  simplified  through  the  application  of  circuit  theory  and  the  representation  of  the 
enclosure  as  a  lossy  inductor.  Inductive  coupling  between  the  enclosure  and  a  stationary  loop 
antenna  placed  coaxially  within  it  was  considered  but  found  to  be  insignificant. 

The  computed  magnetic  field  strength  at  the  center  of  the  enclosure  was  found  to  be  9  dBA/m 
below  the  field  measured  by  the  loop  antenna  at  100  Hz.  The  difference  between  the  computed  and 
measured  results  may  stem  from  various  approximations  made  in  the  model.  Some  key 
simplifications  made  in  the  model  are  as  follows: 

(a)  The  enclosure  was  approximated  as  a  tube. 

(b)  The  loop  antenna  axis  was  placed  coincident  with  the  axis  of  the  enclosure. 

(c)  The  perturbation  of  the  earth's  dc  magnetic  field  in  the  vicinity  of  the  enclosure  due  to  the 
presence  of  the  magnetic  core  was  ignored. 

(d)  The  perturbation  of  the  ac  magnetic  field  of  the  vibrating  enclosure  due  to  the  presence  of  the 
magnetic  core  was  neglected. 

(e)  End  effects  associated  with  both  the  enclosure  and  magnetic  core  were  ignored. 

The  inner  radius  of  the  tubular  enclosure  was  determined  by  equating  its  inside  cross-sectional 
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area  with  that  of  the  actual  enclosure.  This  was  done  since  the  enclosure  was  modeled  as  an 
inductor,  making  it  important  to  equate  the  magnetic  fluxes  passing  through  each  enclosure. 

Nonetheless,  this  analysis  shows  that  the  vibrating  enclosure  produces  a  significant  magnetic 
field.  It  has  been  assumed  that  the  earth's  magnetic  field  lies  in  the  plane  of  vibration  of  the 
enclosure.  Thus  the  computed  magnetic  field  is  a  maximum. 

An  alternative  analysis  can  be  performed  by  solving  Maxwell’s  equations  in  the  rest  frame  of 
the  moving  medium  (enclosure)  through  the  use  of  the  Maxwell-Lorentz  Transformation  [10].  In 
this  method,  the  induced  current  distribution  in  the  enclosure  and  resulting  magnetic  field  are 
determined  in  the  rest  frame  of  the  enclosure  followed  by  a  transformation  of  the  magnetic  field  to 
the  rest  frame  of  the  earth.  This  approach  was  considered  but  appeared  quite  unwieldy. 
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a)  4-in.  Thick  Welded  Aluminum  Enclosure 

b)  Rigid  Antenna  Carriage 

c)  Loop  Antenna 

d)  Vibration  Isolation  System 
(e)  5-in.  Bolted  Base 


antenna  system. 
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Figure  2.  Types  of  vibratory  motion. 

(a)  Translational  motion  in  the  radial  direction. 

(b)  Translational  motion  in  the  axial  direction. 

(c)  Rotational  motion. 
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